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q . ABSTRACT 

!— i ' We report on a search for the optical counterparts of 175 ^tm - selected sources from 

\ the Far-Infrared Background (FIRBACK) survey in the European Large Area ISO 

■ Survey (ELAIS) N2 field. Applying a likelihood ratio technique to optical catalogues 

from the Isaac Newton Telescope - Wide Field Survey (INT-WFS), we found optical 

• i— i , identifications for 33 out of 55 FIRBACK sources in this field. These were then re- 

■ assessed in the light of associations with the ELAIS final catalogue for the the N2 field, 
5-j ' to yield a final set of 31 associations. We have investigated the nature of this pop- 
ulation through a comparison of their observed spectral energy distributions (SEDs) 
with predictions from radiative transfer models which simulate the emission from both 
cirrus and starburst components. We find the far-infrared sources to be 80 per cent 
star bursting galaxies with their starburst component at a high optical depth. The 
resulting SEDs were used to estimate far-infrared luminosities, star formation rates, 
dust temperatures and dust masses. The N2 FIRBACK population is found to consist 
of four suspected ultra-luminous infrared galaxies (ULIRGs) with Lfir ~ 10 12 L 
and SFRfir > 100 MQyr -1 , a number of luminous infrared galaxies (LIRGs) with 
moderate star formation rates and Lfir. ~ 10 L,-:, and a population of low redshift 
quiescently star forming galaxies. We also discuss the implications of these results for 
current evolutionary models. 

Key words: galaxies: evoluntion - galaxies: fundamental parameters - galaxies: 
starburst - infrared: galaxies - surveys 



1 INTRODUCTION 

Evidence available in the 1960s and 1970s from regions 
of active star formation in nearby galaxies, and our own 
Milky Way, suggested that star forming galaxies would be 
strong sources of thermal infrared radiation, and implied 
that regions of dust and molecular clouds should be found 



within them iSunvaev et al.lll978l) . It was then predicted 
that this dust reprocessed redshifted starlight from distant 
galax ies, would form a far-infrar ed background (FIRB) sig- 
nal iPartridee and Peebledll967lh 

With the advent of the Infrared Astronomical Satel- 
lite (IRAS) in the 1980s hundreds of previously undetected 
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galaxies were discovered, which emit u p to 95 per cent o f 
their total luminosity in the infrared dSoifer et al.l I198tI) . 
This indicated significant dust reprocessing of their starlight 
and implied that star formation rates previously calcu- 
lated from rest frame optical/uv luminosities were only 
a lower limit. Follow up redshift surveys, such as QDOT 
llLawrence et al.lll999T) . showed that these galaxi es were far 
more numerous in the past (e.g. lOliver et alll992l) , requiring 
stronger cosmological evolution than had previously been 
predicted by conventional models for the passive evolution 
of galaxies and providing support was building for models 
in which starbursts ar e important evolutionary mechanisms 
jLonsdale et alJll990l) . 

The IRAS results allowed far tighter limits to be 
put on the level of the pr edicted far-infrared background. 
lHacking and Soiferl dl99ll) predicted that the background 
should be detectable by the Cosmic Background Explorer 
satellite (COBE) at 100 /im and noted that discovering its 
exact level would allow further refinement of evolutionary 
models. No-evolution models placed a lower limit on the 
predicted background level, while detection of a strong back- 
ground would indicate the existence of a population of ob- 
jects with significant cosmological evolut ion. The FIRB was 
finally discovered by iPueet et alJ <ll996t) in the whole sky 
survey taken by FIRAS on COBE just as lHacking and Soiferl 
J1991TI had predi c ted. 

IPueet et alJ il99rJl found a signal which could be ex- 
plained by the thermal emission from particles with tem- 
peratures of 20 and 12 K with emissivities of 3 x 10 -6 and 
5 K with an emissivity of 3 x 10 _J . This detection has fu- 
eled a new generation of galaxy evolution models but it is 
clear that to further differentiate between them we must un- 
derstand the nature of the objects which contribute to the 
far-infrared background. 

1.1 The Far— Infrared Background 

The FIRBACK survey was carried ou t jointly with the 
ELAIS team at 175 fan using ISOPHOT jLemke et alll996ft 
on t he Infrared Space Observatory (ISO) llKessler et alJ 
1996). This survey aimed to resolve the FIRB into discrete 
sources and also allow studies of FIRB fluctuations. FIR- 
BACK covered 3 main fields [South Marano, Nl and N2] 
covering a total area of 4 deg 2 . The final source catalogue 
contains 106 sources with fluxes between 180 mjy and 2.4 Jy 
and a supplementar y catalogue containing 90 sources with 
fluxes 135-180 mjy dPole et alJll99Sl . 12001^ . The northern 
fields lie within the previously defined ELAIS Nl and N2 
survey regions whose selection was based primarily on their 
low Galactic cirrus emission (as quantified by IRAS 100 fim 
emission) and high visibility, given ISO's orbit. This choice 
of fields allows multiwavelength follow-up thus providing a 
global view of galaxy evolution. 

There have been several suggestions as to the nature of 
these 175 fim ISO sources. Two s ources were examined i n de- 
tail in the optical and near-IR bv lChaoman et all <|2002T) who 
found two nearby [z < 1) galaxies with characteristic tem- 
peratures of ~30 K and ~50 K and bolometric luminosities 
(40-200 /im) of 10 11 -10 12 L a . These galaxies were classified 
as ULIRGs with morphologies suggesting the early stages of 
mergers adding weight to arguments for interactions spark- 
ing periods of increased star formation. A sample of ISO 
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Figure 1. Distribution of P ra n values for chosen source counter- 
parts. 

175 /im sources were also examined at sub-mm and near-IR 
wavelengths by ISaiina et alJ i2003T) who found a bimodal 
galaxy population: one population of normal star forming 
galaxies at z ~ and a second more luminous population 
at z ~ 0.4 — 0.9. There has been spectroscopic follow-up 
of the brighte st 175 sources in the FIRBACK South 
Marano field JPatris et alj f2003) which were found to be 
nearby (z < 0.3), heavily extincted (A v ~ 3 with extreme 
cases up to A v =7), star forming galaxies with moderate 
star formation rates (a few 10 M0yr -1 ). These were classi- 
fied at LIRGs, with L m ~ 10 11 L . 

We have gone one step further by seeking associations 
for all 55 N2 175^m ISO sources in the FIRBACK sur- 
vey with optical sources from the INT-WFS catalogues 
in four bands (g',r',i',Z), the other ELAIS survey wave- 
lengths (J, H, K, 15, 6.7, 90 urn 20 cm) and IRAS wa ve- 
lengths (60, 100 /im). (Also see iDennefeld et alJlloO^ for 
investigations into sources in the Nl field). In doing so we 
are able to gain a multiwavelength view of the nature of 
these interesting objects. We have compared their SEDs 
with those predicted by the radiative transfer models of 
Efstath iou and Rowan-Robinsonl J2003T) to provide further 
insight into their properties. 

The layout of this paper is as follows: Section [5] de- 
scribes the multiwavelength data, Section [3] discusses the 
methods of associating the multiwavelength data with the 
175 fim ISO sources, provides a summary of the results 
of this process and briefly describes the optical properties 
of the sources. The radiative transfer models are described 
in Section [I] and Secton |H| describes the comparison of the 
sources' SEDs with predictions from the models. The result- 
ing star formation rates, far-infrared luminosities and dust 
temperature estimations are given in Sections |S| and |7| Sec- 
tion |H| provides a summary of current evolutionary models, 
while conclusions and a discussion are provided in Section 
|^| An Appendix provides details of individual associations. 



2 COMPLEMENTARY DATA 

Optica l associations were sou ght in r' data from the INT- 
WFS iMcMahon et alJ 1200 ll) which used the Wide Field 
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Figure 2. Normalised magnitude distribution of the N2 175 
jim ISO source optical counterparts (r — band) (solid line), dashed 
line gives the r'-band magnitude distribution of the INT-WFS 
catalogue for comparison. 
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Figure 3. Redshift distribution for the 30 175 /im sources with 
spectroscopic or photometric redshifts as listed in Table 1. 



Camera (WFC) on the 2.5 m Isaac Newton Telescope at 
the Observatorio del Roque de Los Muchachos (La Palma, 
Spain). The WFC is formed by four 4k x 2k CCDs. The ar- 
rays have 13.5 fim pixels corresponding to 0.33 arcsec/pixel 
at the telescope prime focus and each one covers an area on 
the sky of 22.8 arcmin x 11.4 arcmin. The total sky cover- 
age per exposure for the array is therefore 0.29 deg 2 . Gaps 
between detectors are typically 20 arcsec and optical ob- 
servaions are carried out allowing for a 10 per cent overlap 
between adjacent pointings for photometric purposes. The 
ELAIS-N2 optical field is approximately 9 deg 2 . The survey 
has a typical completeness limit of 23 mag (Vega) in wave- 
bands U, g', r', i' and Z. However, there are no J7-band 
data for our sources the N2 field. All optical magnitudes 
were calculated directly from the INT-WFS data. 

Associations were then sought with sources in the mul- 
tiwavelength ELAIS band-merged catalogue at 6.7, 15 and 
90 nm with associated data at U,g',r',i',Z,J,H,K and 20 cm 
iRowan- Robinson et all I2004 1 ) which by then had become 



available. This was a survey carried out using ISO covering 
a total of 12 deg 2 over five main fields, three in the north 
(Nl, N2, N3) and two in the south (S1,S2). These areas were 
also surveyed at 20 cm with the Very Large Array (VLA) 
and AT. The separate wavebands making up the ELAIS sur- 
vey each comprise an independent survey, the band-merged 
catalogue simply brings together the results of each of their 
final anayses, as described in Section[3] 

Where detections were not found at 15 and 90 /im up- 
per limits were ext r acted directly from th e survey maps 
llVaccari et all 12004 iHeraudeau et alJ 120041) to a 3cr level 
using aperture photometry with an aperture size of 6 arc- 
sec and 90 arcsec respectively. An aperture correction was 
applied to the 15 /im upper limits since 40 per cent of the 
point spread function (PSF) lies outside the aperture. 

IRAS 100 and 60 /im fluxes or upper limits (3a) were 
obtained at the galaxy optical positions using the IRAS Scan 
Processing and Inegration (SCANPI) internet facility 1 to 
provide extra constraints on the source SEDs. 

Associations were also sought in the Submillime- 
tre Common Use r Bolometer Array (SCUBA) 850 fim 
JScott et al.ll2002l) and Max -Planck Millimeter Bolometer 
array (MAMBO) 1200 /im dCreve et alJ 120041) catalogues 
within 10 and 3 arcsec of the optical position of the 175 /im 
counterpart, respectively. However, none were found for any 
of our 175 fim sources. 



IDENTIFICATIONS: 
RESULTS 



METHOD AND 



The positional accuracy of the ISO 175 sources w as de- 
termined using simulated maps bv lDole et alJ ll200ll) . Since 
93 per cent of artificial sources brighter than 180 mjy in the 
simulated maps were recovered within a radius of 50 arc- 
sec, associations were sought for the N2 sources within a 90 
arcsec radius. This ensured the capture within the source 
error circle of the true source position for even faint source 
IDs. Associations were m ade using a l ikelihood ratio tech- 
nique described in full in lMann et alJ <l2002l) . In brief, the 
likelihood ratio is the ratio of the probability of finding the 
true counterpart to a particular source at the source posi- 
tion with that magnitude to the probability of finding an 
object of the same magnitude there by chance. In terms of 
flux this takes the form: 



LR = 



q( f ) e ( x »y) 

n(f) 



(1) 



where q(f) is the flux probability distribution function for 
the source counterparts, e(x,y) is the probability distribu- 
tion of positional offsets between the source and object and 
n(f) is the surface density of objects per unit flux interval. 
e(x,y) was taken to be a gaussian distribution with a — 20 
arcsec. This choice prevents too many objects being assigned 
moderate likelihood ratios, without excluding objects which 
should be considered or including objects which should not. 
The function q(f ) is unknown for a new source popul ation, as 
is the case here, and, as shown in lMann et alJ <fl997l) . an em- 
pirical estimate is very noisy for a small number of sources, 
q(f) was taken to be a constant. The probability, P ra n, that 

1 http:/ /irsa. ipac.caltech.edu/applications/IRAS/Scanpi/ 
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Figure 4. Examples of the affects of the radiative transfer model parameters on the resulting SEDs, see Section [I] Solid black line: 
Ay = 1, fsB — 1, T« = 100 (comparison model). Black dotted line: Ay = 3, /ss = 1, r« = 100. Dash-dot-dot-dot line: t m = 0. Grey 
dashed line: Ay = 1, fsB = 5, t v = 100. Dashed line: ip = 10. Dash-dot line: r = 8 X 10 10 . Long dashed line: tsb = 80 Myrs. Solid grey 
line: t v = 200. 



a fictitious, randomly-placed source would have a likeliest 
assocation with an optical object in the catalogue produc- 
ing a likelihood ratio at least as high as the source-object 
pair in question was then calculated. Thus good associations 
have low values of P ran , corresponding to a high probability 
that the association is correct. The object with the lowest 
Pran was chosen as the source optical ID, with a cut off for 
non-association of P ran = 0.15. This is an arbitrary value 
chosen by insp ection of the P ra n distribution, see Figure 
It is expected jMann et al.ll20o3) that true IDs will cluster 
at low P ran values with a tail of objects reaching to high 
Pran values due to sources without good associations in the 
optical catalogue, as seen in the plot. 

Counterparts were also sough t in the band-merged 
ELA IS multiwavelength catalogue jRowan- Robinson et alJ 
2004) using a search radius of 60 arcsec consistent, with the 
band merging procedure used in the collation of this cat- 
alogue. The individual wavelength ELAIS catalogues were 
merged, taking into account the different positional accura- 
cies of each catalog ue, sequentially as fol l ows. The 15 /xm and 
20 cm catalogues JVaccari et alJ 12004 ICilieei et all Il999l) 
were first i dentified with optical objects f rom the INT-WFS 
catalogue iGonzalez-Solares et alJ 120031) and then merged 
according to their optical positions using a search radius of 
2 arcsec. These were then matched with the 6.7 /im cata- 
logue using a search radius of 5 arcsec and those 6.7 /im 
sources not merged were associated with optical counter- 
parts. These were then merged with the 90 am an d 175 
(im sources jHeraudeau et al.ll2004l : iDole et alJl200 ll) using 
search radii of 30 arcsec and 60 arcsec respectively. Where 
these were matched with more than one catalogue entry the 
less probable assocations were flagged in the catalogue. After 
this those sources that were still unmatched were associated 
with optical counterparts. These optical counterparts were 



sought bv lRowan-Robinson et all i2004l) independently from 
our work and the resulting optical IDs were then compared. 

In summary, we sought associations on the basis of ap- 
plying the likelihood ratio method to the optical data alone, 
then, with the availablility of the final band-merged ELAIS 
catalogue, we reassessed them. Our final optical and ELAIS 
wavelength ID list was drawn up using the following criteria: 

(i) If our optical ID, selected using the likelihood ratio 
method, agreed with the ELAIS optical ID to which the 175 
/im flux had been assigned then both the ELAIS optical ID 
and band-merged data were accepted. 

(ii) If the ELAIS optical ID had a P ran < 0.15 (cut off 
value see earlier) for association with the 175 /jm ISO source 
then the ELAIS optical ID and band-merged data were ac- 
cepted. 

(iii) If our optical ID and ELAIS optical ID were in fact 
the same galaxy but had different entries in the INT-WFS 
catalogue due to the way the image analyser splits up bright 
galaxies into multiple sources, then the band-merged ELAIS 
data and optical ID were accepted. 

(iv) If our optical ID had P ran < 0.15 but no optical 
ELAIS catalogue match then the source has optical and 175 
fj,m data only (it may also have IRAS detections/upper lim- 
its and 60 /im and 90 /mi upper limits) 

(v) If our optical ID had P ra n > 0.15 the source has no 
association 

For sources in the main N2 FIRBACK catalogue (S > 
180 mjy), which contains 27 detections, we made 22 confi- 
dent optical associations with galaxies, 7 of which agree ex- 
actly with the ELAIS optical ID (category i.) and 6 for which 
we accepted the ELAIS optical ID (categories ii. or iii.). Nine 
confident optical associations were made for sources in the 
supplementary catalogue (5* < 180 mjy, 28 sources), 2 of 
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Figure 5. la starburst model parameter confidence region for 
fsB & n d Ay for FN2 008. Possible models parameter combi- 
nations are shown with small crosses (note f$B = 1 — 100 
and Ay = — 3 were allowed), the best-fitting model is 
marked with a square and the la confidence region is shaded. 
A high resolution version of this figure is available at 
http:/ /www. roe. ac.uk/~elt/N2paper Jmages.html . 
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Figure 6. la starburst model parameter confidence region for 
fsB & n d T v for FN2 008. Possible models parameter combi- 
nations are shown with small crosses (note f$B = 1 — 100 
and t v = 50 — 200 were allowed), the best-fitting model is 
marked with a square and the la confidence region is shaded. 
A high resolution version of this figure is available at 
http:/ /www. roe. ac.uk/~elt/N2paper Jmages.html . 



which agreed with the ELAIS optical ID (category i.) and 
2 of which where the ELAIS optical ID was accepted (cat- 
egories ii. or iii.). There also were two sources in the sup- 
plementary catalogue with stellar IDs, one of which agreed 
with the ELAIS optical ID (category i.) and one of which 
was an accepted ELAIS optical ID (category ii. or iii.). The 
data for the identified sources are given in Table 1. and notes 
on individual sources are given in the Appendix, including 
Pran values of any other plausible identifications. 

The identified sources were morphologically classified 
by visual inspection. We used i'-band images from the 
INT-WFS, and follow ed the same classification method as 
iPostman et alJ i2005tl . The galaxies were classified using the 
classical Hubble sequence: E, E/S0, SO, SO/a, Sa, Sa/b, Sb, 
Sb/c, Sc, Irr. For the Advanced Camera for Surveys (ACS) 
on HST, for galaxies with redshift of around unity, visual 
classifications of this kind have a typical random error of 
25 per cent, this reduces to 6 per cent when galaxies are 
classified in two broad cathegories: 1. early-type, 2. spirals 
and irregulars. These errors were estimated from the average 
scatter found in classifications by four ACS team members. 
The majority of our sources are spiral galaxies, we also have 
a number of SO galaxies, 2 sets of interacting galaxies, 2 pairs 
and 2 irregulars, the morphology of 6 sources was unable to 
be deduced from the images. 

Figure [5] shows the optical (r'-band) magnitude distri- 
bution of the N2 175 /im ISO source optical counterparts 
compared to that of the overall INT-WFS catalogue. The 
source IDs are, in general, bright objects with apparent mag- 
nitudes ranging from 21 to as bright as 16. 

Figure |3 shows the redshift distribution of the sources: 
they lie predominantly at z < 0.3 with only two sources 
having z > 0.5. Redshifts are discussed in more detail in 
Section [S] 




Figure 7. icr starburst model parameter confidence region 
for Ay and r v for FN2 008. Possible models parameter com- 
binations are shown with small crosses (note Ay = — 3 
and t v = 50 — 200 were allowed), the best- fitting model is 
marked with a square and the la confidence region is shaded. 
A high resolution version of this figure is available at 
http:/ /www. roe. ac.uk/~elt/N2paper Jmages.html . 



4 RADIATIVE TRANSFER MODELS 

It was initially at tempted to f it the data to the GRASIL 
model SEDs of ISilva et alJ l)l99Sft . These provide six 
templates for nearby galaxies, three starbursts (Arp 
220, M82 and NGC 6090) and three sprials (M100, 
M51 and NGC 6946). However, these templates did 
not describe the data well. Instead radia t ive t rans- 
fer models of [EfstgHiio u and Rowan- Robinson! l|2003F) and 
lEfstathiou et alTi|200nir were employed to generate SED pre- 
dictions for comparison with the data. These are able to 



Table 1. Multi wavelength data for associated sources, fluxes in mjy. 

nb. Sources listed with P ran greater than the threshold 0.15 due to falling into category (iii) (see Section [31 
FN2 004 has P ra n — 1-0 since the associated galaxy is mistakenly classified as a star in the INT WFS catalogue. 

FIRBACK FIRBACK position optical position 175 fim P ran mag 15 fim 6.7 fim 60 fim 90 (tin 100 fim 20 cm z p hot z spec morphological classification 
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Table 2. Parameter values 
starburst 



non-starburst 



Figure 8. \a non-starburst model parameter confidence region 
for Ay and ip for FN2 001. Possible models parameter com- 
binations are shown with small crosses (note Ay = — 10 
and if) = 0.5 — 20.0 were allowed), the best-fitting model is 
marked with a square and the la confidence region is shaded. 
A high resolution version of this figure is available at 
http:/ /www. roe. ac.uk/~elt/N2paper Jmages.html . 



model both cirrus and starburst emission or a combination 
of the two. The models used are briefly described below. 

Ultraviolet to submillimetre emission from stars embed- 
ded in interstellar dust in galaxies known as the 'infrared 
cirrus component' are modelled. The model consists of the 
input stellar radiation field, an interstellar dust model and 
the radiative transfer treatment of the interaction between 
the two which generates the output SED. 

The input stellar radiation field is taken from the 
Galaxy I sochrone Synthesis Spectr al Evolution Library 
(GISSEL. IrBruzual and CharlotJll993l ). This gives the radia- 
tion from a mass of stars from the ultraviolet to near-infrared 
as a function of time. 

Stars form in molecular clouds and at the early stages 
of their lives remain inside them during which time their 
radiation is absorbed by the dust and reprocessed to the 
infrared . During this pha s e of t heir life the model uses the 
code of lEfstathiou et all fcOOd) to compute the radiation 
emitted by stars and dust. The cloud is assumed to disperse 
7.2 x 10 7 years after star formation. Before this, however, 
due to non-spherical evolution of the cloud, a fraction of 
the starlight / may have been able to escape without dust 
absorption. This occurs a time t m after star formation. 

The radiation field within the galaxy is modelled to be 
due to a large number of randomly oriented molecular clouds 
whose average emission, for stars in the age range t m to 
7.2 x 10 7 years, is approximately 

(i - f)st + fs; 

were i s the emission from a giant spherical molecu- 
lar cloud (lEfstathiou et all EoOO) and St is the em ission 
from the stellar population feruzual and Chariot! Il993f) . 
The emission of stars younger than t m is simply S„ and 
those older than 7.2 x 10 7 years St- The fraction of gen- 
eral starlight that escapes the galaxy unattenuated is also 
parameterised in a 'leak' variable. 

The dust temperature is determined by the intensity of 



r 


8 Gyrs 


1 Myr 


tSB 


galaxy age-t m 




galaxy age 


(universe age at z 


- 0.5 Gyrs) 


1~SB 


40 Myrs 




leak 


0.0 


0.1 


f 


1.0 


1.0 


tm 


20 Myrs 


3 Myrs 


A v 


0-3 


0-10 


Tv 


50-200 




1> 


5 


0.5-20.0 


fsB 


1-100 





the stellar radiation field and is characterised in the model 
by the ratio of the bolometric intensity of the radiation field 
to th at in the stellar rad iation field in the solar neighbour- 
hood (Ma this et al.ll983Tl , ip. The star formation rate, (f>,{t), 
is assumed to have an exponential form with a time scale, 
r, and a Salpeter IMF from 0.1 - 125 M . 
— t 

<fi*(t) oc e T 

In the case where a starburst is included the star formation 
rate is modified to 

4>* B {t) (X j>,{t) + fsBC 'Vb 3 ' 

where fsB is the ratio of the star formation rate at the peak 
of the starburst to that at time 0, tsB is the age of the 
galaxy at the start of the starburst and tsb is the exponen- 
tial timescale of the starburst. 

The interstellar dust model of 

ISiebenmorgen and Kruegell l)l992f) is used, which mod- 
els the effects of both small grains and PAHs. The opacity 
of the dust is characterised by Ay, which controls the 
amount of UV to near-IR light absorbed by dust and 
re-emitted in the IR and sub-mm and therefore the ratio 
of the luminosity in these two bands. The model also has 
the ability to increase the dust extinction of the starburst 
component of the galaxy separately via the parameter r„, 
values of 50, 100, 150, 200 are suppported. In the starburst 
case a grid of models was produced over fsB, Av and 
t v and for non-starburst galaxies a grid of models was 
produced over Av and ip. The values of the parameters 
used are given in Tabled 

The effects of each of the model parameters on the 
resulting SED is shown in Figure 0] The solid black line 
shows an SED produced using the fixed parameters (ip — 5, 
r = 8 x 10 9 , t m = 20 Myrs, r SB = 40 Myrs) and f SB = 1, 
Av = 1, and r„ = 100 for comparison. A model with in- 
creased dust extinction of Ay = 3 is shown by the black 
dotted line: this increases the emission at longer wavelengths 
and decreases that in the optical, as would be expected. The 
dash-dot-dot-dot line shows a model with t m = which dra- 
matically changes the shape of the far-infrared peak due to 
the reduction in the amount of dust reprocessing of radiation 
from newly formed stars. However, any significant amount of 
time spent within a molecular cloud produces a far-infrared 
peak of shape similar to that of the other models shown. 
The peak of the SED is also increased significantly with 
an increase in the strength of the starburst, a model with 
fsB — 5 is shown by the grey dashed line. The other param- 
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eters do not affect the resulting SED vastly, hence, there are 
degeneracies and a spread of parameter values may result 
in a very similar SED. This is evidenced by models with 
%l> = 10 (dashed line), r = 8 x 10 10 (dash-dot), t S b = 80 
Myrs (long dashes) and t v = 200 (solid grey line) all hav- 
ing very similar shaped SEDs. The parameter values chosen 
were in the mid-range o f physically viable values as found by 
lEfstathiou et all i200(J) . The r v r ange covers values kno wn 
to fit normal starbursts (t v = 50) ( Efst athiou et al.l200(t) to 
heavily obscured starb ursts such as Arp 220 (t v — 200), this 
range was also used bv lFarrah et al] (I2003T) to fit a number 
of ULIRGs. 

After fitting the models to the multiwavelength data, 
the parameter values for models within the la % confi- 
dence limit of the best-fitting model were examined. The la 
limit for starburst model fitted SEDs typically encompasses 
t v ± 50, Av ± 0.5 and fsB ± 1. There is also a correlation 
between /sb and Av such that models with higher fsB val- 
ues require low Ay values to fit the data. This is due to 
an increase in each of the parameters increasing the height 
of the far-infrared peak (as shown in Figure therefore 
when both the parameters have high values the height of 
the peak is over-estimated. We also note that it is difficult 
to constrain t v since we do not have data in the region of 
the spectrum which it most affects. For non-starburst fitted 
SEDs the confidence region includes AviOA and ip — — 20. 
We also note that since the effect of t[> on the shape of the 
model SEDs is not marked, it is difficult to constrain this 
parameter and the la confidence region encompasses the 
whole allowed parameter range. However, since it has little 
effect on the shape of the model SEDs this will not affect the 
property estimates calculated from the models. Examples of 
these confidence regions for both the starburst (FN2 008) 
and non-starburst (FN2 001) galaxies are shown in Figures 

GDtoOD 



5 SPECTRAL ENERGY DISTRIBUTIONS 

We have compared the observed spectral energy distribu- 
tions (SEDs) with predictions from the radiative trans- 
fer models of Ef s tathio u and Rowan- Robinson (2003) and 
Efst athiou et al] (l200CTt described in Section gl 

Spectroscopic redshifts were used in the fitting of the 
model SEDs, where availabl e from the ELAIS catalogue 
]Rowan- Robinson et al] 12004') or elsewhere. In other cases 
photometric redshifts were drawn from t he ELAIS cata- 
logue , and calculated bespoke by ImpZ iBabbedee et al] 
2004) for the remainder. We note that since photometric red- 
shifts are typically accurate to about 10 per cent in (1 + z), 
for the spectroscopic redshifts we have available this im- 
plies an accuracy of around 0.11 to 0.14. Therefore, given 
that the solutions are based on a maximum of four optical 
bands, the agreement between our photometric and spec- 
troscopic redshifts is reasonable, the largest margin of error 
being about 20 per cent in (1 + z). It is also important to 
note that there are no catastrophic outliers which can oc- 
cur due to degeneracies in the solution space of photometric 
codes. These small redshift errors do not affect the fitted 
SEDs significantly and the resulting errors on derived prop- 
erties are much smaller than those incurred by the la error 



Table 3. Parameters for the starburst radiative transfer model 
SEDs fitted to the N2 175 fim ISO source SEDs. 



FIRBACK 


A v 


fSB 


Tv 


FN 2 002 


0.4 


2.0 


150 


FN 2 003 


0.3 


1.0 


200 


FN 2 004 


0.0 


1.0 


200 


FN2 005 


0.0 


2.0 


100 


FN2 007 


0.2 


2.0 


200 


FN 2 008 


0.0 


2.0 


150 


FN 2 010 


0.3 


3.0 


200 


FN2 011 


0.5 


3.0 


200 


FN2 012 


0.0 


3.0 


100 


FN 2 016 


0.3 


2.0 


200 


FN 2 017 


0.1 


1.0 


200 


FN2 018 


0.0 


2.0 


200 


FN2 019 


0.0 


1.0 


200 


FN 2 021 


0.0 


1.0 


150 


FN 2 022 


0.1 


1.0 


150 


FN 2 025 


0.8 


3.0 


150 


CFN2 029 


1.1 


10.0 


200 


CFN2 030 


0.9 


2.0 


200 


CFN2 034 


1.0 


1.0 


50 


CFN2 038 


0.0 


1.0 


150 


CFN2 039 


0.0 


1.0 


100 


CFN2 044 


0.0 


1.0 


150 


CFN2 047 


0.0 


1.0 


50 


CFN2 049 


0.1 


5.0 


50 



Table 4. Parameters for the quiescent radiative transfer model 
SEDs fitted to the N2 175 A*m ISO source SEDs. 



FIRBACK 


A v 


i> 


FN 2 000 


10.0 


20.0 


FN 2 001 


0,1 


17.0 


FN2 015 


1.1 


20.0 


FN 2 020 


3.7 


0.5 


FN 2 023 


0.2 


0.5 


CFN2 036 


1.0 


0.5 



on the fitted SED \ 2 value which are given for all derived 
quantities. 

The N2 175 /xm ISO source SEDs are shown in Figures 
rHand lTni They appear to fall into two population categories, 
the first being star bursting galaxies with their starburst re- 
gion at a high optical depth (fitted model parameters shown 
in Table El and SEDs plotted in Figure an d the second 
dust obscured galaxies with quiescent star formation (fitted 
model parameters shown in Table [I] and SEDs plotted in 
Figure EHJ. 



6 STAR FORMATION RATES AND 
FAR-INFRARED LUMINOSITIES 

Star formation rates for the sour ces were calculated using 
both an optic al [/-band estima te dCram et al.lll993l and a 
FIR estimate ]Kennicuttlll998h . given in equations [5] and [3] 

SFRc7(M > 5 M Q ) = -100 Moyr- 1 (2) 

V ~ ; 1.5 x 10 22 W Hz" 1 ° y V ' 

SFRfir = 4.5 x 10- 44 L F i R (ergs s" 1 ) Meyr" 1 (3) 
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Figure 9. N2 ISO 175 /im sources fitted model SEDs. Best-fitting starburst models (solid line) and cirrus models (dashed line) are 
shown for each source. 
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Figure 9 - continued N2 ISO 175 fim sources fitted model SEDs. Starburst models (solid line) and cirrus models (dashed line) are 
shown. 
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Figure 9 - continued N2 ISO 175 fim sources fitted model SEDs. Starburst models (solid line) and cirrus models (dashed line) are 
shown. 
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Figure 10. N2 ISO 175 fim sources fitted model SEDs. Cirrus models (solid line) and starburst models (dashed line) are shown. 



Where Lu is the [/-band luminosity, estimated through in- 
terpolation of the SED models and Lfir is the integrated 
infrared luminosity over 8-1000 fim. These both assume a 
Salpeter IMF. The chosen IMF affects the estimated star 
formation rates since we integrate over only a portion of the 
luminosity of the galaxy's constituent stars. The proportion 
of stars emitting with luminosities in the range in question 
will depend on the initial mass function, since it describes 
the initial distr ibution of stellar lu minosities over the range, 
as discussed bv lMann et alJ J2002I) . The star formation rates 
for individual sources can be seen in Table |5] The distribu- 
tion of FIR star formation rates and a comparison between 
those estimated using both the FIR and [/-band luminosi- 
ties can be seen in Figure ITT1 There is a large spread in star 
formation rates with a number of galaxies with SEDs sug- 



gestive of a modest starburst. This implies that these galax- 
ies are of low mass since their existing stellar population is 
dominated by the small starburst at infrared wavelengths. 
Using star formation rates estimated for SED models falling 
in the k x 2 range (see Section the [/-band and FIR 
SFRs are estimated to be accurate to 25 per cent and 50 
per cent respectively. It can be seen in the lower pannel of 
Figure 1111 that the FIR star formation rate estimates are 
significantly higher, typically by a factor of 4, than those 
using the U-band estimate. This is to be expected for such 
highly obscured objects. 

The integrated far-infrared luminosities over the wave- 
length range 8-1000 fj,m were estimated. Results can be seen 
in Table and a distribution in Figure WI\ These are esti- 
mated to be accurate to approximately 25 per cent. 
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Table 5. Star formation rates and FIR luminosities of the star- 
burst galaxies (first block) and cirrus galaxies (second block). 



FIRBACK 


t/-band SFR 


FIR SFR 


Lfir 




z 




Meyr" 1 


M yr- 1 


L 






FN2 002 


36.9 


93.1 


4.87 x 


10 11 


0. 


17 


FN 2 003 


3.8 


3.6 


1.83 x 


10 10 


0. 


03 


FN 2 004 


9.4 


5.4 


2.85 x 


10 10 


0. 


03 


FN2 005 


131.4 


169.1 


9.18 x 


10 11 


0. 


26 


FN2 007 


22.3 


39.5 


2.07 x 


10 11 


0. 


12 


FN2 008 


24.4 


27.3 


1.46 x 


10 11 


0. 


11 


FN2 010 


21.7 


88.9 


4.68 x 


10 11 


0. 


17 


FN2 011 


104.8 


391.9 


2.04 x 


10 12 


0. 


51 


FN2 012 


6.3 


13.1 


7.07 x 


10 lu 


0. 


11 


FN2 016 


12.6 


25.6 


1.33 x 


10 11 


0. 


13 


FN2 017 


32.5 


28.1 


1.41 x 


10 11 


0. 


18 


FN2 018 


8.4 


12.2 


6.53 x 


10 10 





15 


FN2 019 


1.7 


1.0 


5.22 x 


10 9 


0. 


.03 


FN2 021 


7.9 


4.8 


2.59 x 


10 10 


0. 


10 


FN2 022 


32.6 


22.3 


1.18 x 


10 11 


0. 


.17 


FN2 025 


5.8 


39.7 


2.07 x 


10 11 





20 


CFN2 029 


0.1 


7.2 


3.79 x 


10 10 


0. 


15 


CFN2 030 


3.2 


16.8 


8.60 x 


10 10 


0. 


15 


CFN2 034 


138.2 


296.6 


1.48 x 


10 12 


0. 


55 


CFN2 038 


25.6 


11.6 


6.19 x 


10 10 





.14 


CFN2 039 


7.7 


4.1 


2.23 x 


10 10 


0. 


07 


CFN2 044 


0.4 


l.S 


9.41 x 


10 9 


0. 


05 


CFN2 047 


63.2 


27.2 


1.49 x 


10 11 





25 


CFN2 049 


49.1 


252.8 


1.39 x 


10 12 





51 


FN2 000 


2.8 


7.9 


3.97 x 


10 10 


0. 


.03 


FN2 001 


9.3 


1.7 


8.37 x 


10 9 





03 


FN2 015 


29.1 


18.8 


9.58 x 


10 10 


0. 


17 


FN2 020 


3.0 


7.5 


2.56 x 


10 10 


0. 


15 


FN2 023 


1.1 


0.1 


3.87 x 


10 8 


0. 


02 


CFN2 036 


4.3 


2.5 


8.82 x 


10 9 





10 



Table 6. Greybody parameters and dust masses for the starburst 
galaxies (first block) and cirrus galaxies (second block). 



FIRBACK T (K) (3 M dust (M Q ) 



FN2 002 


40 


1. 


.4 


8. 





X 


10 7 


FN 2 003 


36 


1. 


.5 


4. 


.3 


X 


10 6 


FN2 004 


41 


1 


.9 


1. 


.4 


X 


10 6 


FN2 005 


40 


1 


.9 


2. 


.9 


x 


10 7 


FN2 007 


42 


1. 


.5 


2. 


.5 


x 


10 7 


FN2 008 


39 


1 


.9 


6. 


.5 


x 


10 e 


FN2 010 


45 


1 


.4 


5. 


.4 


x 


10 7 


FN2 Oil 


41 


1 


.4 


2. 


.5 


x 


10 8 


FN2 012 


39 


1. 


.9 


3. 


.0 


X 


10 e 


FN2 016 


12 


1 


. i 


2. 


.1 


x 


10' 


FN2 017 


36 


1 


,1 


3. 


.9 


X 


10 7 


FN2 018 


39 


2 





2. 


.4 


X 


10 6 


FN2 019 


41 


1. 


.9 


2. 


.5 


X 


10 5 


FN2 021 


41 


1 


.9 


9. 


.9 


X 


10 r> 


FN2 022 


41 


1 


.6 


1. 





X 


10 7 


FN2 025 


40 


1 


.3 


4. 


.5 


X 


10 7 


CFN2 029 


19 


1 


.2 


5. 


6 


X 


10 f ' 


CFN2 030 


38 


1. 


.3 


2. 


.5 


X 


10 7 


CFN2 034 


35 


1 


.1 


6. 


.5 


X 


10 8 


CFN2 038 


40 


1 


.9 


2. 


6 


X 


10 6 


CFN2 039 


40 


1. 


.9 


9. 


.2 


X 


10 s 


CFN2 044 


41 


1 


.9 


i. 


2 


X 


10 6 


CFN2 047 


41 


1 


.7 


6. 


6 


X 


10 e 


CFN2 049 


44 


1 


.5 


6. 


6 


X 


10 7 


FN2 000 


31 


1 


.9 


7. 


.2 


X 


10 e 


FN2 001 


30 


1 


.9 


1. 


7 


X 


10 e 


FN2 015 


31 


1 


.9 


1. 


.4 


X 


10 7 


FN2 020 


17 


1. 


.8 


1. 


6 


X 


10 8 


FN2 023 


16 


2 


.0 


2. 


.0 


X 


10 e 


CFN2 036 


17 


1 


.8 


5. 


1 


X 
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7 DUST TEMPERATURES AND MASSES 



The dust in the galaxy radiative transfer models is at a 
range of temperatures, however, we computed an emission- 
weighted estimate via the fitting of a grey body to the far- 
infrared portion of the SED (from the peak of the SED to 
1000 (im) as follows:. 

F v oc — (4) 

evr — 1 

In equation 2] T is the dust temperature which controls the 
position of the maximum of the long wavelength SED, j3 is 
the frequency dependence of the dust grain emissivity and 
controls the gradient of the long wavelength slope of the grey 
body curve. Fitting over this wavelength range allows us to 
constrain both the temperature and (3, without reaching the 
far-infrared where the dust may go optically thick and no 
longer be well described by a greybody model. The derived 
values for the sources can be seen in Table |S] Using the 
method described above (Section^ for the starburst models 
these are estimated to be accurate to T±3 K and /3±0.5, in 
the non-starburst case these estimates are accurate to T±14 
K and /3±0.1. 

The distribution of the dust temperatures for the N2 
175 /jm ISO source galaxies is shown in Figure 1131 The 
starburst gala xy dust temperature s are slightly higher than 
those found bv lsaiina et~af] J2003) (~ 20 - 30 K) for a sam- 



ple of Nl I SO 175 urn sources but have similar emissivity 
coefficients: ISaiina et al.l <l2003l) found /3 ~ 1.5 — 1.7. 

Dust m asses were estima ted using the submillimetre 
prescription jHildebrandlll983l) at a rest wavelength of 450 
/im: 



Mdust 



(1 + z)n(v r )B(v r ,T dust ) 



= ao67 (2^Toir)' 



(5) 
(6) 



where vq and v r are the observed and rest frame frequen- 
cies respectively, S VQ is the observed flux at z/o, estimated 
through interpolation of the SED models, and B{v r ,Td U st) 
is the Planck function in the rest frame. Tdust and (3 are 
the dust temperature and dust grain emissivity estimated 
via the greybody fitting. «(i/ r ) is known as the mass ab- 
sorption coefficient and has units of m kg -1 . This approach 
assumes that the galaxy is optically thin at submillimetre 
wavelengths. Resulting dust masses, calculated at a wave- 
length of 450 /im (accurate to approximately 25 per cent), 
are shown in Table HJ We note that the dust masses follow 
the trend of the far-infrared luminosities and star forma- 
tion rates, those with high SFRs and Lfir having larger 
dust masses. 
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Figure 13. Distribution of dust temperatures for the ISO 175 
lira sources, see Table IS1 
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Figure 11. Top window:Far— infrared star formation rate distri- 
bution (MQyr- 1 ) for the ISO 175 fim sources, see Table 151 Bot- 
tom window: Far-infrared starformation rate estimates vs. u-band 
starformation rate estimates. 
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Figure 12. Distribution of far— infrared luminosities (Lq) for the 
ISO 175 /im sources, see Table 151 



8 EVOLUTIONARY MODELS 

The detection of the far-infrared background in 1996 led to a 
re-evaluation of galaxy evolution mo dels. One of the e arliest 
new models to emerge was that of iBurigana et all (I1997T) 
which used the FIRB to place constraints on models for the 
cosmic star formation history and therefore the production 
of the background signal. This paper concluded that stellar 
nucleosynthesis is able to account for a FIRB intensity not 
far from the upper limit set by the COBE detection. It also 
suggested substantially higher star formation rates at higher 
redshifts than suggested by optical data. 

It is now fairly well established that infrared-selected 



galaxies undergo extremely high rates of evolution compared 
to measurements at other wavelengths over the approximate 
redshift range < z < 1. Chary and Elbaz (2001) models 
imply that about 80 per cent of the cosmic infrared back- 
ground was produced between < z < 1.5 compared to 
only 30 per cent of the 850 (j,m background. Also the dust 
enshrouded star formation rate peaks at z ~ 0.8 at a value 
of ~ 0.25 M yr _1 Mpc" 3 and at least 70 per cent of this 
star formation takes place in infrared luminous galaxies with 
Lir > 10 11 hp,. In fact ~ 7 5 per cent of 15 fim sources exam- 
ined bv lElbaz et alJ i2002l) were found to be in this category 
of luminous infrared galaxies and have star formation rates 
of ~ 100 M©yr -1 . The comoving density of infrared light 
due to these galaxies was over 40 times larger at z ~ 1 than 
it is today. 

Models of the co smic infrared background produced by 
iLaeache et alJ (2003'). that are able to fit number counts and 
unresolved anisotropies imply that the luminosity function 
must change dramatically with redshift. These also show 
rapid evolution of high luminosity sources up to z=l. It has 
been suggested that this strong evolution is due to bimodal 
star formation, one long lived quiescent, passive phase and 
one enhance d starburst phase triggered by interactions an d 
mergers (e.g.lLagache et aljEooa T lFranceschini et al]l200ll) . 
iFranceschini et alTl20o"ll) also suggest that the interpreta- 
tion of evolution with redshift could be biased due to an 
increase in the probability of detecting a galaxy during an 
active phase caused by both an increase of the rate of in- 
teractions, in part due to the geometry of the expanding 
universe, and an increase of the infrared luminosity due to 
more abunda nt fuel available i n the past. 

However. lChapman et all l)2002|) point out that at z < 1 
the coldest and most dusty galaxies will have the great- 
est flux densities for an equivalent infrared luminosity. At 
z > 1 the FIRBACK 175 fim detection will lie bluewards 
of the peak in the restframe black body spectrum and dust 
temperatures of < 100 K will no longer have significant ef- 
fects on the observed flux. Therefore at higher redshifts, the 
flux limited FIRBACK survey will be biased towards cooler 
galaxies. 

Models such as IFranceschini et all (1200 ll) and 
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ICharv and Elbazl |2001) do not explicitly incorporate 
such cold luminous galaxies whic h are preferenti al ly se- 
lected in our survey for z < 1. IChapman et alJ (2002) 
examined two of the 175 /jm selected ISO sources and 
indentified them as unusually cold ULIRGs with dust 
temperatures of ~ 30 K which is significantly lower than 
the usual ULIRG dust temperature of ~ 50 K. Cold, dusty 
yet luminous objects suggest large masses of dust heated 
at a moderate intensity. The existence of such sources may 
seriously affect the interpretation of SCUBA sources if they 
are numerous at high-z. Since many models rapidly evolve 
a hot ULIRG population to the local luminosity function 
they may predict a false redshift distribution if these cold 
sources are not considered. 



9 DISCUSSION AND CONCLUSIONS 

We have found bright galaxy optical identifications for 31 
out of the 55 N2 ISO 175 /im sources and stellar identifica- 
tions for 2 sources. 17 of those with galaxy optical IDs were 
associated with sources from the band-mergerd ELAIS cat- 
alogue. We then compard the spectral energy distributions 
of the source s with predictions from the radi a tive trans- 
fer models of [E fstathio u and Rowan- Robinson! (2003) and 
lEfstathiou et alJ (|200(tl ~ 

22 of the 175 fim sources have not been confidently asso- 
ciated with optical counterparts and therefore their nature 
is still unknown. However, only 4 of these sources are in the 
main FIRBACK catalogue, the remaining 18 are listed in 
the complementary catalogue which contains detections at 
lower signal to noise and therefore they may just be spuri- 
ous detections. However, these could also be due to sources 
having multiple IDs or them being associated with heavily 
extincted objects which were not in the optical catalogue. 

24 galaxies have been classed as star bursting. They 
have a range of far-infrared luminosities from 10 9 Lq to 
10 12 Lq corresponding to a range of star formation rates 
from ~ 0.1 to ~ 400 M,@yr~ . Although those at the lower 
end of this range are not forming stars at a dramatic rate 
their small starburst component is dominating over the ex- 
isting population at far-infrared wavelengths. This implies 
that the 175 fim FIRBACK sample contains a significant 
number of low mass galaxies. The majority of these galaxies 
are spirals and a number are SO galaxies. 

Sources FN2 005, 011, CFN2 034 and CFN2 049 have 
Lfir ~ 10 12 Lq and high star formation rates, we therefore 
suggest that these are in the low redshift tail of the rapidly 
evolving ULIR G population. T his is in agreement with the 
findings of lSaiina et al.l <2003h who found one sixth of their 
175 jUm sample to fall into this category. We also suggest that 
those galaxies with moderate star formation rates and far- 
infrared luminosities of ~ 10 11 Lq at moderate z (> 0.12) 
are LIRGs (FN2 002, 007, 008, 010, 016, 017, 022, 025 and 
CFN2 047). Those with z < 0.12 are more normal galaxies 
perhaps with lower masses (FN2 003, 004, 012, 018, 019, 
021, CFN2 029, 030, 038, 039, 044). 

The estimated dust temperatures of the starburst galax- 
ies are 30-49 K which is relatively low compared to that of 
standard ULIRGs which have Td ~ 50 K. The most lumi- 
nous of these galaxies also have large dust masses adding 
weight to the suggestion that cold, dusty, luminous objects 



have large masses of dust heated at a moderate intensity 
JChapman et alJl2002l) . 

A further 6 sources have non-starburst SEDs. These 
sources have far-infrared luminosities of 10 8 — 10 10 Lq and 
star formation rates of < 30 Moyr -1 . We suggest that these 
are quiescently star forming galaxies (FN2 000, 001, 015, 
020, 023, 036). These galaxies have low dust temperatures of 
16—31 K and dust masses of 10 6 -10 8 M Q . 3 of these galaxies 
are spirals, 1 is an intereacting galaxy pair, 1 is either an 
irregular or a spiral with tidal tail and the morphology of 
the remaining sour c e was unable to be deduced. 

iLaeache et alJ (2003) predicted 62 per cent of 175 /im 
sources would have z < 0.25 and the remainder 0.8 < z < 
1.2. We find only one 175 /im source with S > 180 mjy with 
z > 0.25 and therefore have S> 62 per cent with z < 0.25. 
However, our mix of higher redshift ULIRGs and lower red- 
shift starburst galaxies is in good qualitative agreement 
with the predictions of the [ La g ache e^iL| i2003 f) starburst 
templ ates as adopted bv lSaiinaet^ITT2OT3T ~ lLagache et alJ 
(2003) also state that to gain agreement between their mod- 
els and observations they require a cold local population 
which is what has been found. 

There has also been recent investigation into the sub- 
millimetre properties of the local univ erse through the 
SCUBA Local Universe Galaxy Survey iDunne and Ealesl 
1200 ll iDunne et al]|200ch . This survey found a population 
of galaxies with Lfir = 10 10 — 10 11 Lq and dust masses 
of 10 7 - 10 s M . The population was fitted with a two- 
component dust model to discover cold dust temperatures 
of ~ 21 K and hot dust temperatures of 30 — 50 K. It is clear 
that the local 175 ^m FIRBACK population has very similar 
properties to these galaxies. The sources of the far-infrared 
background are therefore offering a stepping stone between 
high redshift SCUBA galaxies and the local universe. 
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11 APPENDIX 

Notes on the multiwavelength associations for individual 
sources. Images of each of the sources are shown in Figure 

m 

FN2 - 000: Identification with a bright optical galaxy 
(Pran = 0.026, r' magnitude = 18.5, z spec — 0.03, 
irregular or edge-on spiral with tidal tail). This is 
also an ELAIS source (ELAISC15_J163734. 4+405208) op- 
tical ID with detected emission at 15, 6.7, 90 fim 
and radio detected emission at 20 cm. There are 
two other ELAIS sources (ELAISC15_J163729 . 3+405248, 
ELAISC15_J163731. 3+405156) within the FIRBACK error 
ellipse with optical IDs with associated P ran values of 0.604 
and 0.246. 

FN2 - 001: Identification with a bright optical galaxy 
(Pran = 0.002, r' magnitude = 18.3, z phot — 
0.10, edge-on Sb/c). This is also an ELAIS source 
(ELAISC15_J163507. 7+405929) optical ID with detected 
emission at 15, 6.7, 90 pm and radio detected emission at 
20 cm. This also has an IRAS detection at 60/im. 

FN2 - 002: Identification with an optical galaxy (P ra n = 
0.078, r' magnitude = 19.3). An ELAIS source ID 
(ELAISC15_J163608. 1+410507) has an optical ID with the 
same galaxy, however, they are different entries in the 
INT catalogue caused by the splitting up of bright ob- 
jects by the image analyser. The ELAIS ID has P ran = 
0.086, z spec ~ 0.17, it is an irregulare galaxy with a 
tidal tail and detected has emission at 15, 6.7 and 90 
pirn. This also has an IRAS detection at 60 /im. There 
are two other ELAIS sources (ELAISC15_J163608 . 4+410529, 
ELAISC15_J163611. 6+410452) within the FIRBACK error 
ellipse with optical IDs with associated P ra n values of 0.885 
and 0.983. 

FN2 - 003: Identification with a bright optical galaxy 
(Pran = 0.001, r' magnitude = 17.5, z apec = 0.03, Sa). This 
is also an ELAIS source (ELAISC15_J163525 . 1+405542) op- 
tical ID with detected emission at 15, 6.7, 90 /jm and radio 
detected emission at 20 cm. This also has an IRAS detection 
at 60 /jm. 

FN2 - 004: This ID is a small optical galaxy (P ran = 0.079, 
r magnitude = 20.2). However, there is also present an- 
other large bright galaxy (r magnitude = 16.4) which is 
mistakenly classified as a star so is forcably not associated 
and assigned a P ran of 1.0. Previous associations allowing 
ISO 175 fim sources to be associated with stars chose this 
ID which matches with the optical ID of an ELAIS source 
(ELAISC15_J163401. 8+412052) which has z spec = 0.03 and 
detected emission at 15, 6.7, 90 /im and radio emission at 
20 cm and has the morphology of a SBa. This also has an 
IRAS detection at 60 /im. 

FN2 - 005: Identification with bright optical galaxy 
(Pran = 0.011, r magnitude = 18.7, z spec — 0.26, Sa). 
This is also an ELAIS source (ELAISC15_J163242 . 4+410847) 
optical ID with detected emission at 15, 6.7, 90 fim 
and radio emission at 20 cm. This also has an IRAS 
detection at 60 /im. There is another ELAIS source 
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Figure 14. Optical INT images of the N2 ISO 175 /an sources. FIRBACK source position marked with a diamond (axes show 
distance from source position in arcsec), optical ID marked with a cross, ELAIS positions marked with small circles, ELAIS optical IDs 
marked with pluses. A 90 arcsec circle is plotted to show the region in which optical IDs were searched for. Figure available from 
http: / /www. roe. ac.uk/~elt /N2paper Jmages.html 



(ELAISC15_J163244. 6+410911) within the FIRBACK error 
ellipse with an optical ID with an associated P ran of 0.315. 

FN2 - 006: Identification with what appears on the im- 
age to be a pair of interacting galaxies (P ran = 0.001, 
r' magnitude = 18.2, no redshift). An ELAIS source 
(ELAISC15_J163506. 1+411038) appears to be on the same 
object, however, is does not have an optical ID. There is an- 
other ELAIS source (ELAISC15_J163503. 5+411137) within 
the FIRBACK error ellipse. 

FN2 - 007: Identification with bright compact source with 
two nuclei (P ran = 0.022, r magnitude = 18.0, z sp( , c — 0.12). 
This is also an ELAIS source (ELAISC15_ J163546 . 9+403903) 
optical ID with detected emission at 15, 6.7, 90 /im 
and radio emission at 20 cm. This also has an IRAS 
detection at 60 /im. There is another ELAIS source 
(ELAISC15_J163545 . 5+403825) within the FIRBACK error 
ellipse with an optical ID with an associated P ran of 0.998. 

FN2 - 008: Identification with an optical galaxy (P ran = 
0.053, r' magnitude = 18.4). However, the 175 /jm flux 
has been assigned in the ELAIS catalogue to a source 
(ELAISC15_J163548. 0+412829) with a different optical ID. 
This ID has P ran = 0.084 for association with the FIRBACK 
source and is therefore a possible identification. The ELAIS 
source has z spec — 0.14 and detected emission at 15, 6.7, 90 
/im and radio emission at 20 cm, it is an Sb galaxy. This 
also has an IRAS detection at 60 /im. 

FN2 - 009: This source does not have a confident opti- 
cal identification. It is best associated with a galaxy with 
Pran — 0.440, this does not agree with the optical ID for the 
ELAIS source (ELAISC15_J163359. 2+405303, P ran = 0.880, 
r' magnitude=18.8) to which the 175 /im flux has been as- 
signed, this source also has detected emission at 15, 6.7, 
90 /im and radio emission at 20 cm. This source is confi- 
dently identified as a star when stellar associations are al- 
lowed with P ran = 0.142 (the threshold for these associa- 
tions is P ran = 0.4), however, this still does not agree with 
the ELAIS optical ID. 

FN2 - 010: Identification with a pair of interact- 
ing galaxies, an irregular with a tidal tail and an Sa 
(P ran = 0.114, r' magnitude = 18.8). An ELAIS source 
(ELAISR163536+411726) ID has an optical ID with the same 
galaxy, however, they are different entries in the INT cat- 
alogue caused by the splitting up of bright objects by the 
image analyser. The ELAIS source has z apec — 0.17 and de- 
tected emission at 90 /im and radio emission at 20 cm. This 
also has an IRAS detection at 60 /im. 

FN2 - 011: Identification with an optical galaxy (P ran = 
0.118, r' magnitude = 19.8, z p hot = 0.51, Sa?). However, 
the 175 /im flux has been assigned in the ELAIS catalogue 
to a source (ELAISR163809+405839T) with a different galaxy 
ID (r magnitude = 18.6) and detected emission only in the 
far-infrared and radio. This has a P ran value of 0.547 for 
association with the FIRBACK source, it is therefore not 
an acceptable association. 

FN2 - 012: Identification with an optical galaxy (P ran = 
0.006, r' magnitude = 18.5). However, the 175 /jm flux 



has been assigned in the ELAIS catalogue to a source 
(ELAISC15_J163412. 0+405652) with a different galaxy ID 
(r magnitude = 20.0). This ID has P ran = 0.088 for 
association with the FIRBACK source and is therefore 
a possible identification. The ELAIS source has z 3pec = 
0.14 and detected emission at 15, 6.7, 90 /im and ra- 
dio emission at 20 cm, it is possibly an Sa. There are 
two other ELAIS sources (ELAISC15_J163417. 4+405710, 
ELAISC15_J163417. 9+405653) within the FIRBACK error 
ellipse with optical IDs with associated P ran values of 0.990 
and 1.000. 

FN2 - 013: This source does not have a confident optical 
identification. It is best associated with an optical galaxy 
(P ran = 0.260, r' magnitude = 20.4). This is also an ELAIS 
source (ELAISC15_J163406 . 5+405106) optical ID with de- 
tected emission at 15 /im. 

FN2 - 014: This source appears in the ELAIS 90 and 175 
/im unassociated catalogue (ELAIS-FBK175_N2_014). It does 
not have a confident optical identification. 

FN2 - 015: Identification with a pair of interacting galax- 
ies, an Sa/b and an Sa, (P ran = 0.015, r magnitude 
= 18.6, z spec — 0.17). This is also an ELAIS source 
(ELAISC15_J163607. 7+405546) optical ID with detected 
emission at 15, 6.7, 90 /im and radio emission at 20 cm. 

FN2 - 016: Identification with an edge-on op- 
tical galaxy whose morphological type is unclear 
(P ran = 0.119, r' magnitude = 20.0). An ELAIS source 
(ELAISC15_J163423. 9+405410) ID has an optical ID with 
the same galaxy, however, they are different entries in the 
INT catalogue caused by the splitting up of bright objects 
by the image analyser. This source has z 3pec — 0.13 and 
detected emission at 15, 6.7 /im and radio emission at 20 
cm. This also has IRAS detections at 60 and 100 /im. There 
are two other ELAIS sources (ELAISC7_J163421+405413, 
ELAISC7_J163423+405506) within the FIRBACK error 
ellipse one with an optical ID with P ran = 1.000, the other 
has no optical ID. 

FN2 - 017: Identification with a pair of interacting galax- 
ies, an Sa and an edge-on spiral, (P ran = 0.030, r' mag- 
nitude = 18.5, Zpi lot = 0.18). The 175 /im flux is as- 
signed to an ELAIS source (ELAISR163442+410759) with 
a galaxy ID with associated P ran = 1.000, it is very faint 
and far from the source. There is, however, another ELAIS 
source (ELAISR163445+410817) present, with detected emis- 
sion only in the radio, in the error ellipse without an optical 
ID which is only about 5arcsec from our FIRBACK optical 
ID. 

FN2 - 018: Identification with an optical galaxy (P ran = 
0.037, r' magnitude = 19.5, z phot = 0.15, SO/a). How- 
ever, the 175 /im flux has been assigned in the ELAIS 
catalogue to a source (ELAISC15_J163334 . 1+410139) with 
a different galaxy ID with P ran = 0.998 (r magnitude = 
20.5), it is therefore not a possible ID. This also has an 
IRAS detection at 60 /im. There is another ELAIS source 
(ELAISR163333+410112) within the FIRBACK error ellipse 
whose galaxy ID has a P ran = 1.000 
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FN2 - 019: Identification with an optical galaxy (P ran = 
0.001, r' magnitude = 18.1, z phot = 0.10, Sb). The ELAIS 
source (ELAISR163716+404825) to which the 175 /im flux 
has been assigned does not have an optical ID, it has de- 
tected emission at 6.7/im and 20 cm. This also has an IRAS 
detection at 60 /im. 

FN2 - 020: Identification with an optical galaxy (P ran = 
0.086, r'magnitude = 20.2, z phot = 0.15). The ELAIS source 
(ELAISC15_J163242 . 7+410627) to which the 175 /im flux has 
been assigned is fainter (r magnitude = 22.1) and further 
from the source than our ID and therefore has a higher P ran 
value of 0.990. 

FN2 - 021: Identification with an optical galaxy (P ran = 
0.001, r'magnitude = 17.6, z p hot = 0.10). However, the flux 
of this FIRBACK source is below the ELAIS catalogue 5a 
limit of 223 mjy and therefore as it has no associations at 
other ELAIS wavelengths it is omitted from the catalogue. 

FN2 - 022: Identification with an optical galaxy 
(P ran = 0.019, r' magnitude = 19.1). The ELAIS source 
(ELAISC15_ J163708 . 1+412856) to which the 175 fim flux has 
been assigned has an optical identification with a bright op- 
tical galaxy with r' magnitude = 18.8, it has P ran = 0.140 
and is therefore a possible identification for the FIRBACK 
source. This ELAIS source has z spec = 0.17, is an Sa and 
detected emission at 15, 6.7 fim and radio emission at 20 cm. 
There is another ELAIS source (ELAISC7_J163709+412832) 
within the FIRBACK error elipse, however, its optical ID is 
stellar and therefore it has P ran = 1.0. 

FN2 - 023: Identification with an optical galaxy (P ran = 
0.018, r' magnitude = 19.0, z phot = 0.02, Sb/c). However, 
the flux of this FIRBACK source is below the ELAIS cata- 
logue 5a limit of 223 mjy and therefore as it has no asso- 
ciations at other ELAIS wavelengths it is omitted from the 
catalogue. 

FN2 - 024: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 
the ELAIS catalogue 5a limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

FN2 - 025: Identification with an optical galaxy (P ran = 
0.119, r' magnitude = 20.2, z phot = 0.20, Sa).The ELAIS 
source (ELAISC7_J163628+404757) to which the 175 /mi 
flux has been assigned has optical ID classed as a galaxy, 
however, it looks stellar and has a P ran of 0.974 so 
is not a possible ID. There is another ELAIS source 
(ELAISC7_J163633+404749) within the FIRBACK error 
elipse, however, its optical ID is stellar and therefore it has 
P ran = 1.0. This also has IRAS detections at 60 and 100 
/im. 

FN2 - 026: This source does not have a confident op- 
tical ID. It is associated with a galaxy with r' magni- 
tude 22.4 with a P ran = 0.906. There is an ELAIS source 
(ELAISC15_J163615. 7+404759) within the FIRBACK error 
ellipse, however, it has a stellar optical ID and therefore has 
P ran = 1.0. 

CFN2 - 027: This source does not have a confident optical 
ID. It is associated with a galaxy with r' magnitude 20.7 
with a Pran — 0.305. This does not agree with the ELAIS 
source (ELAISR163703+412425) optical ID which is a galaxy 
with P ran = 0.999. 

CFN2 - 028: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 



the ELAIS catalogue 5a limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

CFN2 - 029: Identification with an optical galaxy 
(P ran = 0.070, r' magnitude = 19.6, z phot = 0.15). 
The 175 /im flux has been assigned to a fainter galaxy 
(ELAISR163418+410729, r' magnitude = 22.2) which is 
further from the FIRBACK position than our ID, it 
therefore has a high P ran of 1.0. This also has an IRAS 
detection at 60 /im. There are three other ELAIS sources 
(ELAISR163419+410641, ELAISC15_J163421 .4+410622, 

ELAISR163422+410648) within the FIRBACK error ellipse, 
however, none of them have optical IDs. 

CFN2 - 030: Identification with an optical galaxy (P ran = 
0.086, r' magnitude = 20.3, z p hot — 0.15, SO/a). However, 
the flux of this FIRBACK source is below the ELAIS cata- 
logue 5a limit of 223 mjy and therefore as it has no asso- 
ciations at other ELAIS wavelengths it is omitted from the 
catalogue. However, it has IRAS detections at 60 and 100 
/im. 

CFN2 - 031: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 
the ELAIS catalogue 5a limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

CFN2 - 032: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 
the ELAIS catalogue 5a limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

CFN2 - 033: This source does not have a confident optical 
ID. It is associated with a galaxy with r' magnitude 21.8 
with a P ran = 0.934. 

CFN2 - 034: Identification with an optical galaxy (P ran = 
0.030, r' magnitude = 19.4, z phot = 0.55, SO/a). However, 
the flux of this FIRBACK source is below the ELAIS cata- 
logue 5a limit of 223 mjy and therefore as it has no asso- 
ciations at other ELAIS wavelengths it is omitted from the 
catalogue. 

CFN2 - 035: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 
the ELAIS catalogue 5a limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

CFN2 - 036: Identification with an optical galaxy (P ran = 
0.119, r'magnitude = 19.4, z p hot = 0.10, Sa). However, the 
flux of this FIRBACK source is below the ELAIS catalogue 
5a limit of 223 mjy and therefore as it has no associations at 
other ELAIS wavelengths it is omitted from the catalogue. 
However, it has IRAS detections at 60 and 100 /im. 

CFN2 - 037: This source does not have a confident optical 
ID. It is associated with a galaxy with r' magnitude 22.1 
with a P ran = 0.812. When associations are allowed with all 
objects this source is identified as a star with P ran = 0.183. 
However, this still does not agree with the stellar optical ID 
of the ELAIS source (ELAISR163812+405452) assigned the 
175 fim flux. 

CFN2 - 038: Identification with an optical galaxy 
(P ran = 0.023, r' magnitude = 19.3, z spec — 0.14, 
Sb/c). This is also the optical ID for an ELAIS source 
(ELAISC15_ J163431 . 5+412246) with detected emission at 15 
fim and has an IRAS detection at 100 /im. 
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CFN2 - 039: Identification with an optical galaxy 
(P ran = 0.001, r' magnitude = 17.2, z apec = 0.07, 
Sb). This is also the optical ID for an ELAIS source 
(ELAISC15_J163613. 6+404230) with detected emission at 
15, 6.7, 90 fim and radio emission at 20 cm. This also has 
an IRAS detection at 60 /im. 

CFN2 - 040: This source does not have a confident 
optical ID. It is best associated with a galaxy with r' 
magnitude 20.5 and P ran = 0.192. The ELAIS source 
(ELAISC15_J163641. 1+413131) to which the 175 /im flux 
has been assigned has a galaxy optical identification 
with P ran = 0.562. There is another ELAIS source 
(ELAISC15_J163648. 1+413134) within the FIRBACK error 
ellipse with a galaxy optical ID with P ran = 0.393. 

CFN2 - 041: This source does not have a confi- 
dent optical ID, however, its best association agrees 
with the optical identification of the ELAIS source 
(ELAISC15_J163433. 6+405953) to which the 175 fim 
flux has been assigned. This galaxy has r' magni- 
tude = 19.2 and P ran = 0.311, the ELAIS source 
has detected emission at 15 /im and 20 cm. There 
are two other ELAIS sources (ELAISR163427+405936T, 
ELAISC15_J163430. 1+410055) within the FIRBACK error 
ellipse, one has a galaxy optical ID with P ran = 0.998 and 
the other does not have an optical identification. 

CFN2 - 042: This source does not have a confident galaxy 
optical association. The flux of this FIRBACK source is also 
below the ELAIS catalogue 5<j limit of 223 mjy and there- 
fore as it has no associations at other ELAIS wavelengths it 
is omitted from the catalogue. 

CFN2 - 043: This source does not have a confident op- 
tical identification when associations are only allowed with 
galaxies (P ran = 0.437). However, when all objects are pos- 
sible associations it is matched with a star (P ran = 0.056) 
which is also the ELAIS source (ELAISC7_J163546+404929) 
ID to which the 175 /im flux has been assigned. This ELAIS 
source has detected emission at 6.7 /im. 

CFN2 - 044: Identification with an optical galaxy (P ran = 
0.053, r' magnitude = 20.0, z phot = 0.05, SO/a). The ELAIS 
source (ELAISC15_J163730. 4+404542) to which the 175 /im 
flux has been assigned has a galaxy ID with a P ran value 
of 0.995. This also has an IRAS detection at 60 /im. There 
is another ELAIS source (ELAISR163728+404533) within the 
FIRBACK error ellipse, however, this source does not have 
an optical ID. 

CFN2 - 045: This source does not have a confident galaxy 
optical association. The flux of this FIRBACK source is also 
below the ELAIS catalogue 5cr limit of 223 mjy and there- 
fore as it has no associations at other ELAIS wavelengths it 
is omitted from the catalogue. 

CFN2 - 046: This source does not have a confident optical 
association. The flux of this FIRBACK source is also below 
the ELAIS catalogue 5<j limit of 223 mjy and therefore as it 
has no associations at other ELAIS wavelengths it is omitted 
from the catalogue. 

CFN2 - 047: Identification with an optical galaxy (P ran = 
0.070, r' magnitude = 19.0). The 175 /im flux has been as- 
signed to an ELAIS source (ELAISC15_J163449. 5+412048) 
with a galaxy optical ID with P ran — 0.125 which is there- 
fore a possible association. This source has z spec = 0.25, is 
an SO/a and detections at 15, 6.7 /tm and 20 cm. There is an- 
other ELAIS source (ELAISC15_J163451 . 9+411944) within 



the FIRBACK error ellipse with a galaxy optical ID with 
P ran = 0.911. 

CFN2 - 048: This source does not have a confident optical 
identification, is is best associated with a galaxy with P ran = 
0.920. The ELAIS source (ELAISC15_J163739. 2+405643) to 
which the 175 /tm flux has been assigned does not have an 
optical identification. 

CFN2 - 049: Identification with an optical galaxy (P ran = 
0.092, r' magnitude = 20.3, z phot = 0.18). An ELAIS source 
(ELAISC15_J163741. 3+411913) ID has an optical ID with 
the same galaxy, however, they are different entries in the 
INT catalogue caused by the splitting up of bright objects by 
the image analyser. The ELAIS source has detected emission 
at 15 and 90 /im and is possible an Sa. This also has an 
IRAS detection at 60 /im. There is another ELAIS source 
(ELAISC15_J163738. 9+41 1840) within the FIRBACK error 
ellipse with a galaxy optical ID with P ran = 0.911. 

CFN2 - 050: This source falls between two WFS chips. 

CFN2 - 051: This source does not have a confident optical 
identification, it is best associated with a galaxy with P ran = 
0.866. The ELAIS source (ELAISC15_J163616. 2+41 1426) to 
which the 175 /im flux has been assigned has a stellar optical 
ID and therefore is not a possible ID for the FIRBACK 
source. 

CFN2 - 052: This source does not have a confident optical 
identification, it is best associated with a galaxy with P ran = 
0.466. The ELAIS source (ELAISC15_J163403 . 0+410350) to 
which the 175 /im flux has been assigned has a galaxy op- 
tical ID with P ran = 0.962. There are two other ELAIS 
sources (ELAISC7_J163407+410219, ELAISR163410+410326) 
within the FIRBACK error ellipse one with a stellar ID and 
one with a galaxy ID, however, both have P ran = 1.0. 

CFN2 - 053: This source does not have a confident optical 
identification, it is best associated with a galaxy with P ran = 
0.441. The ELAIS source (ELAISR163621+412244) to which 
the 175 /u,m flux has been assigned has a galaxy optical ID 
with P ran = 0.989. 

CFN2 - 054: This source does not have a confident optical 
identification. When associated with all optical objects it is 
confidently associated with a star with P ran = 0.215, the 
ELAIS source (ELAISC7_J163658+411417) assigned the 175 
/tm flux also has a stellar ID with P ran = 0.350 which is 
below the limit of 0.4 and therefore a possible association. 



